Hamiltonian tomography in an access-limited setting without state initialization 
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We propose a scheme for the determination of the coupling parameters in a chain of interacting spins. This 
requires only time-resolved measurements over a single particle, simple data post-processing and no state initial- 
ization or prior knowledge of the state of the chain. The protocol fits well into the context of quantum-dynamics 
characterization and is efficient even when the spin-chain is affected by general dissipative and dephasing chan- 
nels. We illustrate the performance of the scheme by analyzing explicit examples and discuss possible exten- 
sions. 
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It is often the case in the dynamics of interacting many- 
body systems that a specific and desired effect is achieved by 
means of an appropriately designed set of coupling strengths. 
Such effects are frequently very sensitive to even small de- 
viations from the required pattern of interaction parameters, 
which may result in a dramatic deterioration of the perfor- 
mances. Moreover, in many cases the couplings are assumed 
to have been pre-engineered by a third party and it would be 
highly desirable, from a practical viewpoint, to test if we have 
been provided with the proper set of parameters before run- 
ning a protocol. In other words, it would be important to have 
a diagnostic and non-invasive routine which allows one to in- 
fer the pattern of interaction strengths in a quantum many- 
body system with a high degree of accuracy. In essence, we 
would require the performance of a "Hamiltonian tomography 
scheme". This can be seen as a variation of quantum process 
tomography which, together with state H] and detector 
tomography yO, allows for the complete characterization of 
quantum dynamics. All of them have found experimental ver- 
ification. Our scheme greatly reduces the resources necessary 
to estimate the coupling parameters of the interaction model. 

Reconstructing the form of a given but undisclosed one- 
and two-spin Hamiltonian has raised the interest of the physics 
community [4] and, very recently, an intriguing proposal has 
been put forward for the TV-particle case (5|]. Within the 
frameworks of these investigations, protocols able to find the 
coefficients characterizing the interaction Hamiltonian have 
been developed. However, the initialization of the state of 
the system is required in each of them. In addition, a com- 
plete set of relevant eigenvalues of the interaction model has 
to be known a priori or should be determined in an adept 
way, which may require the enforcement of strong conserva- 
tion laws on the class of Hamiltonians that can be tested |5|]. 
These requirements are in general difficult to be met or unnec- 
essarily limiting. In this Letter, we use an approach based on 
the "information flux" (IF) [6] to investigate Hamiltonian to- 
mography performed with minimal access to the many-body 
system and without the necessity for initial preparation |0]. 
Moreover, we stress another remarkable advantage in the pro- 
tocol we suggest: the Hamiltonian to study does not need to 
commute with the total spin- excitation number, i.e., we do not 
require that the total number of excitations in the system is 



preserved. In clear contrast even with classical schemes for 
inverse problems in vibration [8], our method does not rely on 
the prior knowledge of a set of eigenvalues of the Hamilto- 
nian. Information about the coupling coefficients is found via 
time-resolved single- spin measurements without state initial- 
ization of the system, which is distinctive and original, com- 
pared to what has previously been done 10, [2]. 

The difference with respect to quantum process tomogra- 
phy is also evident: there, the initialization of the whole sys- 
tem in a set of relevant states and the performance of state 
tomography after the action of the process are required. From 
the reconstructed output density matrices, one can then infer 
the completely positive map corresponding to the process it- 
self. On the other hand, in our scheme, we just need to mea- 
sure a single element of a multipartite register at various times; 
no condition on the state of the rest of the system is imposed. 
The time evolution of the expectation value of operators act- 
ing on that single spin can be extracted from the acquired data 
and the complete set of coupling coefficients of the Hamilto- 
nian can be reconstructed from it. 

To fix the ideas and clearly elucidate the main features of 
our study, we start from a simple excitation-preserving class 
of interaction models. We consider a linear chain of N spin- 
1/2 particles, mutually coupled via the nearest-neighbor anti- 
ferromagnetic XX Hamiltonian ^\ 

N-l 

Hi=J2 UXiXi+i + Y& +1 ). (1) 

i=l 

Here, Ji > is the interaction strength between spins i and 
i + 1 while Xi, Yi and Z{ denote the x, y and z-Pauli ma- 
trix of spin i, respectively. While it is important to remark, 
at this stage, that our method can be adapted to a larger class 
of Hamiltonians, as discussed later on, we clarify here that the 
choice of Hi is made simply to provide an immediate intuition 
of the protocol through a significant example. The dynamics 
encompassed by Hi will be analyzed via the IF fl6|]. In partic- 
ular, within the framework of explicitly limited accessability 
stated above, we focus our attention on the evolution of qubit 
1. We thus move to the Heisenberg picture and consider the 
dynamics of the Pauli operators of this spin under the action 
of Hi . From now on, time-evolved operators will be indicated 
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FIG. 1 : (a) Sketch of the scheme for Hamiltonian tomography with- 
out state initialization, where M is the measurement performed on 
spin 1 and other spins are not accessible; (b) Simulated dynam- 
ics of (Xi(t)), sampled at steps of Jt = 7r/25, under the action 
of Td with N = 6 and {Ji/J} = {1.02, 1.26, 0.94, 1.36, 0.72} 
(case 1), {Ji/J} = {1.49,0.80,1.02,0.69,1.28} (case 2), and 
{Ji/J} = {1.30, 0.80, 1.23, 0.75, 0.96} (case 3). The correspond- 
ing fits are performed using the trial function in Eq. ([3]), which is in 
excellent agreement with the behavior of the data. 

as G(t) = U^OU with = e-W™ 1 *. A straigthforward 
calculation, based on the use of the operator expansion theo- 
rem, leads to the following decompositions of X(t) and y(t) 
over sets of TV- spin operators 111011 

X 1 (t)= ai(*)Xi+ Z x \a 2 (t)%+ • • +a N (t)Z 2 ■ -R N ], 

with Rn = Xn and Sn = Yn (Rn = Yn and Sn = Xn) 
for odd (even) N. The time-dependent parameters a^(£)'s and 
f3i(tys(i = 1 , . . . , N) are functions of the coupling strength set 
{ Ji}. For instance, we have a\{t) = X^o[(^)V^]^ w * m 
the recurrence formula 5® = (— I) 3 '[Jj-\5^Z-P + Jofif+i^]* 
Jo = Jn = and the initial conditions 6^ = (1) for j ^ 1 
(j = 1). Therefore, a.\(t) depends on the full set { Ji}. 

Focusing, for the sake of argument, on the first of Eqs. ©, 
if we initialize spin 1 in an eigenstate of Xi, \± x )i = (|0) ± 
|l))i/\/2, we get (Xi(t)) = ±a 1 (t). Analogously, if the ini- 
tial state of spin 1 is \± v ) 1 = (|0> ±i|l»i/V2, we would 
obtain (Yi(i)) = ±/3i(t). It is easy to see that, for the case of 
Hi , the recurrence formulas that determine on (t) 's and (t) 's 
are exactly the same, so that ot\(t) = Pi(t). This is due to 
the symmetric role played by the XiX i+ i and YiYi+i terms 
in Eq. (OQ). As stated above, a\(t) depends on the full set 
{ Ji}. Therefore, in order to obtain information about all the 
coupling strengths within Wi, we need to determine the func- 
tional behavior of the expectation value of a single one-qubit 
operator, such as {X\ (t)), and make explicit the relation con- 
necting it to Ji's. As we show later, this can be done via a 
simple post-processing step of the measured data. 

At first sight, it might seem that spin 1 should be properly 
prepared. However, this is clearly not the case. In fact, we just 
need to sample (Xi(t)) at successive instants of time, so that 
each X -projection we perform on spin 1 prepares it into the 
desired class of states. From that point on, (Xi(t)) = ±ai(t) 
holds rigorously. The necessity of iterated state initializa- 
tions is thus bypassed. A sketch of the scheme is presented 



in Fig.QJa). A second important observation is that, by hav- 
ing decoupled the evolution of (cr = X, Y) from the 
explicit influences of the expectation value of operators in- 
volving spins from 2 to N, the initial state of the rest of the 
system might be completely arbitrary and unknown. To the 
best of our knowledge, this feature is unique to our method. 
As no information is required on the dynamical aspects of the 
rest of the chain, our Hamiltonian tomography is performed 
with only minimal invasiveness on the many-body system. 

In order to show the efficiency of the method and clarify its 
working principles, it is worth addressing a few explicit exam- 
ples. We have generated random sets of coupling parameters 
(for a chain of N = 6 spins) taken from a uniform distri- 
bution in the range [0.5J, 1.5 J], where J is an arbitrary con- 
stant. We have then simulated X -measurements on spin 1 in 
a way so as to get a 25-point sample of {X\ (t)) for each case, 
taking t G [0, tt/J] 111 ill at steps of 7r/(25J). The elements 
of three of such samples are shown in Fig. \Hb). As a post- 
processing stage of our analysis, we now need to fit the points 
within each sample with a proper functional form which we 
take as a linear combination of trigonometric cosine functions 
of unknown amplitudes and frequencies. The choice of such 
basis of functions is not arbitrary and is somehow induced by 
the interference nature of the mechanism behind information- 
propagation across a spin-chain, as discussed in 111 211 . More- 
over, this form can also be inferred from the functional form 
of ai(t) in the particular case of sets {Ji} allowing perfect 
state transfer 111311 . For the case at hand, we find that the trial 
function 

a i Ty) CO = A cos(u A t) + B cos(u; B t) + C cos(c^) (3) 

is in excellent agreement with the behavior of the simulated 
data, as shown in Fig. [T] By equating the amplitudes A,B,C 
and frequencies uja,b,c to the functions of J^'s entering into 
ai(t), we have estimated Ji/ J (i = 1, .., 5) to be within 0.1% 
of the values listed in the caption of Fig. [TJb) for each of the 
cases shown. It is remarkable that, differently from previous 
proposals QSQ, the energy spectrum of the coupling Hamil- 
tonian is required at no stage of the protocol. For an ideal 
unitary evolution, less points within each sample are actually 
sufficient to estimate the parameters. For instance, we have 
obtained J^'s with a good precision by considering only 10 
points and light computational effort for the fit. However, the 
plots presented here include 25 points, as this helps in opti- 
mizing their visualization. Although the analysis above pro- 
vides clear evidence that the accuracy of the protocol is almost 
insensitive to the particular choice of J^'s, we have explicitly 
checked this feature by simulating the performance of the to- 
mography protocol for several randomly generated sets, eval- 
uating the average error associated with the retrieval of the 
corresponding coupling parameters. We have studied chains 
of up to N — 8 spins and found an average relative error al- 
ways smaller than 0.3% Ill4ll . 

As previously mentioned, the method can be extended to a 
more general Hamiltonian model. In fact, let us consider the 
following interaction model, which does not preserve the total 
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number of excitations in the system 

N-l 



n 2 =Y. (Jx,iXiX i+1 + Jy^y 1+1 ). 



(4) 



The evolution of X\ and Y\ under H 2 can still be written as 
in Eqs. 0. The coefficients a^'s and /3k s depend, in this 
case, on two disjoint and alternate sets of parameters J a /s 
(this result has been recently exploited in Ref. 111511 ). For in- 
stance, a/e's (/3k s) depend only on Jx,fc's with even (odd) k 
and Jy,/c's with odd (even) k. Actually, the same recurrence 
formulas used above also hold in the present case. We can 
thus perform the tomographic protocol twice: first we con- 
sider |d= £C ) 1 as the initial state of the first qubit and evalu- 
ate (Xi(t)). Then, we estimate (Yi(t)) from the initial state 
\± y ) v In this way we obtain information on both sets of pa- 
rameters and we can reconstruct the complete set {J a ,i}- It 
is important to stress that this excitation-non-preserving case 
cannot be analyzed by means of protocols such as the ones in 
Ref. [5], which critically rely on the condition of excitation- 
conservation. 

The tomographic scheme we propose is able not only to 
provide information on the unitary dynamics of the elements 
of the chain but also to estimate the effects of incoherent cou- 
pling of the system with an environment. In what follows, 
we describe how the influence of dissipation and dephasing 
on a spin-chain can be retrieved from our formal apparatus. 
For the sake of argument, we consider again model Hi al- 
though, we remark, the results will be valid for the more gen- 
eral Hamiltonian H 2 as well. We assume weak-coupling con- 
ditions between each spin of the chain and its own bath, mod- 
eled as an ensemble of bosonic modes. We consider the ef- 
fects of both amplitude and phase damping channels over the 
chain. Using the operator-sum representation, the evolution 
of an initial state p c of the whole chain under the effect of one 
of such channels is given by q c (t) — ^ K^(r)p c K^ (r). 

Here, {K^(r)} is the set of time-dependent Kraus operators 
such that (t)K^(t) = t and r is the interval dur- 

ing which the channel is acting QjJ] . The formal description of 
a single-spin amplitude damping process in a bath in equilib- 
rium at a finite temperature is described by the set {K- 1 (r ) } = 
{A°,ji if, if}, wh ere if = Vp(|0)(0| + e~ r -/ 2 1 1> <1 1 ), 
A\ = y/pq-e -r^lQXil, 

|1>(1|), if = v^l-pXl-e- 
l)/(2n + 1) and n the average phonon number of the bath, 
assumed to be the same for each spin. For a dephasing chan- 
nel, on the other hand, we have {Kf (r)} = D}} with 
4° = y / (l + e-^)/2i i , D\ = y/(l-e-^ T )/2Zi. T 
and 7 are the rates of amplitude and phase damping, respec- 
tively. Our approach is to intersperse the unitary evolution 
Q c (t) = Up c W (induced by Hi) and the non-unitary dynam- 
ics (resulting from the integration of the Lindblad equation 
corresponding to a given noise channel), each lasting for small 
time interval At. We randomly select the spin upon which ap- 
ply the operator-set {A^} or {D^ }. The results of the sim- 
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FIG. 2: (a) Simulated dynamics of (Xi(t)) in the presence of inco- 
herent coupling of the system with an environment. We have con- 
sidered N = 6, 7/ J = 0.5, r/ J = 0.2, n = 0.01, and the sets 
{Ji/ J} listed in the caption of Fig. [Tib). The corresponding fits are 
performed using the trial function in Eq. ((TJ), which is in excellent 
agreement with the behavior of the data; (b) Simulated dynamics of 
(Xi(t)) in the same setting, with a finite number of measurements 
rimeas = 500 and the set {Ji} corresponding to case 1 of Fig. [Tib). 



ulated measurements are then averaged over the collection 
of noise-occurrence patterns (runs), in order to guarantee the 
faithful unraveling of the open quantum dynamics. Finally. 
Hamiltonian tomography is performed. An elegant and effec- 
tive description of open dynamics can be given in terms of IF 
formalism considering that the action of a set of Kraus opera- 
tors on<jj is obviously given by 0(r) = ^2 K^(t)OK^(t). 
Therefore, for a dephasing channel acting on qubit i, we have 



Xi 



r X r 



(5) 



It is immediate to recognize that this results in the change 
ai(i) -> e-^ T ai(t) [(3i(t) -> e~ 7T "/%(*)] in the decompo- 
sition of X\(t) [yi(t)\ in Eqs. ©, with all the other terms 
unmodified. The description in Eqs. will thus remain for- 
mally the same, together with all the qualitative results pre- 
sented in the unitary case. Computationally, the analysis per- 
formed with this method is faster than the one based on den- 
sity matrix evolution (the computational time grows as TV 2 ). 
For an amplitude damping channel, we have 



Xi = e 2 Xi , y 7 
Zi = (l- 



Nry. 



rr )(2p-l)i, + e- Fr Z,. 



(6) 



The decomposition of X\(t) [j)i(t)] after the action of the 
channel is therefore no more restricted to the operator set used 
in Eqs. (|2]) but involves a larger one. For our numerical study, 
we have taken j/J = 0.5, Y/J = 0.2 and n = 0.01. In 
Fig. 0a) we present the simulated points obtained for a chain 
of 6 qubits with the sets { Ji] previously considered and 100 
runs. We have found that the trial function 

a i nc) CO = [A cos(out) +23 cos(cj^) +C cos(^ c t)]e"^ (7) 

with 7 an effective rate depending on 7 and T, is in excellent 
agreement with the behavior of the simulated data, as shown 
in Fig. Ha). Surprisingly, the amplitudes A,B,C and frequen- 
cies uja,b,c are the same as in the ideal case. The only net 
effect of noise in (X\(t)) is the damping of the oscillations. 
Also in the presence of incoherent coupling of the system with 
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an environment, our Hamiltonian tomography protocol works 
well. We have proved it by obtaining A,B,C and uja,b,c from 
the fits in Fig.Oa) and estimating Ji / J. The results are within 
4% of the original values, for each of the cases analyzed. 

In order to predict the performance of our tomography pro- 
cess in a way so as to be closer to realistic conditions, we have 
considered the error due to the finite number of measurements 
^meas performed to evaluate, each time, the required expecta- 
tion values. We have estimated the parameters Ji 's in the three 
cases above, including noise, for n meas = 500. The results are 
within 9% of the expected values, for each of the cases an- 
alyzed. This error can be reduced by increasing the number 
of measurements per sampling time. A set of simulated out- 
comes of the measurements is presented in Fig. [2b). 

Finally, for the sake of completeness, we have analyzed the 
case with additional (unknown) spurious terms in the coupling 
model, such as local magnetic fields along the z-axis or inter- 
action terms proportional to ZiZi+\. This would get our study 
closer to a true experimental situation where unwanted "engi- 
neering" defects could affect a Hamiltonian. In the limit of 
small influences ( ~ 0.1 J) of the additional terms on Hi or 
H2 and for only a finite number of measurements being per- 
formed, our protocol can estimate the parameters J^'s with a 
10% error. This value is comparable to the one obtained with- 
out spurious couplings and an equally finite sampling, which 
shows that the effects of the additional terms is very small. 

We have proposed a scheme for the tomography of a wide 
class of interaction Hamiltonians. Our method is designed to 
work in a scenario of restricted accessability to the compo- 
nents of a spin chain. It allows the identification of coupling 
parameters through the temporal dynamics of a single spin. 
As no initial state preparation is necessary, measurements can 
be performed by interspersing the system's evolution. Besides 
data acquisition, only a simple post-processing step is neces- 
sary: no conservation law associated to the interaction or a 
priori knowledge on the state of the system is required. Even 
when a spin chain is affected by environmental influences, our 
Hamiltonian tomography remains possible and reliable. In or- 
der to widen the class of Hamiltonians that can be assessed 
with our scheme, including the case of external local fields ap- 
plied to the system, we can speculate various strategies, such 
as dynamical changes of measurement basis. A general the- 
ory for Hamiltonian tomography will be an exciting extension 



of the work presented in this Letter. Given the crucial role 
that proper coupling patterns have in the interference effects 
behind quantum many-body phenomena, non-demanding di- 
agnostic methods are important tools which need to be devel- 
oped. Our proposal contributes to this task in a significant 
way. 
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